In this paper, sound propagation in shallow water environments is studied. A Quadratic Finite Element Method (QFEM) by means of functional technique in conjunction with narrow angle parabolic model of hyperbolic wave equation is applied by solving wave equation in depth dependent and range independent fluid-solid media. The capability of FEM based method for modeling of complex solid-fluid boundaries in bottom and also the random surface is employed for formulating of depth operators in Parabolic Equation (PE) method. After simplifying of some integrals in QFEM, the governing system of equations are obtained. Due to existence of mass matrix in the final system, implicit methods seem to be reliable tools for solving such algebraic system of equations. Hence, Crank-Nicolson approach is used for estimation of Solution. Checking of proposed method efficiency is determined by some standard test problems. The obtained results show acceptable agreement with the physical behavior of wave propagation nature but for finding of more accurate and capable methods, it requires to develop this method for wide angle PE in associated with more numerical tricks especially in bottom interfaces.